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Abstract 

We  show  the  dielectrophoretic  actuation  of  single-crystal  diamond  nanomechanical  devices 
using  gradient  radio-frequency  electromagnetic  forces.  Both  cantilever  and  doubly  clamped 
beams,  fabricated  using  our  angled-etching  fabrication  technique, 1  are  demonstrated,  with  op¬ 
eration  frequencies  ranging  from  a  few  MHz  to  ~50  MHz.  Frequency  tuning  and  parametric 
actuation  are  also  studied. 

Owing  to  its  large  Young’  modulus,  excellent  thermal  properties,  small  thermo-elastic  dissi¬ 
pation,  single-crystal  diamond  is  an  ideal  candidate  for  realization  of  high-frequency  (f)  and  high 
quality  factor  ( Q )  nanoscale  mechanical  resonators.  These  devices  are  of  interest  for  realization 
of  stable,  high /  •  Q  product,  RF  oscillators  and  inertial  sensing  applications.2  Recently,  diamond 
nanomechanical  resonators  embedded  with  luminescent  crystalline  defects  (color  centers),  have 
also  been  explored  as  a  promising  platform  for  applications  in  quantum  information  science  and 
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technology  (QIST).3  Of  these,  negatively  charged  nitrogen-vacancy  (NV-)  color  center  is  of  par¬ 
ticular  interest  since  it  can  be  used  as  a  spin-qubit  with  optical  read-out.  Importantly,  NV-  electron 
spin  state  can  have  long  coherence  times  (milliseconds)4  even  at  room  temperature  and  it  can  be 
prepared  and  manipulated  using  microwave  and  radio  frequency  (RF)  fields  to  drive  transitions  be¬ 
tween  these  electron  and  nuclear  spin  sublevels.5-7  Recently,  it  has  been  proposed  that  strain-fields 
can  also  be  used  to  manipulate  the  spin-state  of  NV,8,9  which  resulted  in  renewed  interest  of  QIST 
community  in  single-crystal  diamond  micro-electromechanical  systems  (MEMS). 10-12  For  exam¬ 
ple,  coupling  between  an  NV-  center  and  a  mechanical  resonator  could  enable  high  fidelity  control 
of  NV“  spin  state  via  rapid  adiabatic  passage, 10  and  potentially  the  remote  coupling  of  NV  cen¬ 
ters  via  mechanics.3  Furthermore,  mechanical  resonators  may  enable  coherent  coupling  between 
dissimilar  systems  with  degrees  of  freedom  possessing  dramatically  different  properties  and  en¬ 
ergy  scales.  To  realize  such  systems,  it  is  essential  to  achieve  a  reliable  method  for  controlling 
and  reading  the  dynamic  action  of  monolithic  single-crystal  diamond  mechanical  resonators.  To 
date,  only  a  few  reports  have  demonstrated  mechanical  resonators  fabricated  in  single-crystal  dia¬ 
mond,  13-18  with  recent  results  demonstrating  the  coupling  of  NV-  centers  to  bulk  phonon  modes 10 
and  micron-scale  mechanical  resonators. 11,12  Furthermore,  the  only  true  actuator  realized  in  single¬ 
crystal  diamond  relied  on  heavily  boron  doped  substrates  in  order  to  achieve  electrostatic  actuation 
for  switching  applications.13,14  Such  schemes  are  not  easily  applied  in  QIST  applications  since 
doping  inevitably  compromises  the  diamond  material  properties  and  stability  of  color  centers. 

Dielectrophoretic  actuation 19  has  recently  been  used  to  actuate  and  transduce  motion  of  na¬ 
noelectromechanical  systems  (NEMS)  and  applied  to  achieve  mechanical  resonance  tuning,20  and 
coherent  control  of  classical  mechanical  resonators,21  as  well  as  to  study  cavity  electromechan¬ 
ics,22  and  nonlinear  mechanics.23  In  our  approach,  single-crystal  diamond  cantilevers  (Figure  1(a)) 
and  doubly  clamped  nanobeams  (Figure  1(b))  are  placed  between  metal  electrodes  and  are  driven 
by  fringing  radio-frequency  (RF)  field  (Figure  1(c))  with  frequency  resonant  with  mechanical  res¬ 
onance  of  the  devices.  Our  numerical  modeling  indicates  that  keeping  small  vertical  distance  from 
electrodes  to  nanobeam  is  of  crucial  importance  for  efficient  actuation  (Figure  1(d)). 
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In  contrast  to  capacitive  and  piezo-electric  actuation  approaches,  that  may  require  coating  or 
doping  of  diamond  nanomechanical  structure,  gradient  RF  force  does  not  require  any  modification 
to  the  mechanical  resonator  itself.  Therefore,  it  preserves  diamond’s  unique  material  properties 
(Young’s  modulus,  low  mass  density,  large  thermal  conductivity,  and  low  thermoelastic  dissipa¬ 
tion)  thus  minimizing  dissipations  and  resulting  in  devices  with  large  mechanical  quality  factors. 
Of  particular  interest  for  QIST  applications,  dielectrophoretic  actuation  scheme  maintains  the  in¬ 
tegrity  and  stability  of  NV  centers  embedded  within  diamond  nanomechanical  structures,  that  are 
known  to  be  sensitive  to  fabrication  imperfections  and  surface  terminations.24  Additionally,  we 
note  that  electrodes  employed  for  actuation  may  also  double  for  on-chip  delivery  of  microwave 
power  to  NV-  center.  Therefore,  we  believe  that  dielectrophoretic  actuation  is  an  ideal  choice  for 
investigation  of  NV-mechanical  coupling. 

Our  diamond  cantilevers  and  doubly  clamped  beams  are  fabricated  using  fabrication  technique- 
angled-etching-that  we  have  recently  developed, 1  which  is  shown  in  Figure  2(a).  Briefly,  angled- 
etching  employs  anisotropic  oxygen  plasma  etching  at  an  oblique  angle  to  the  substrate  surface, 
yielding  suspended  triangular  cross-section  nanobeams  directly  from  single-crystal  bulk  diamond 
substrates.  Nanomechanical  resonators  fabricated  by  angled-etching  are  easily  integrated  with 
dielectrophoretic  actuation.  The  distance  from  the  bottom  apex  of  the  suspended  nanobeams  to 
the  substrate  was  carefully  minimized,  since  the  small  vertical  distance  from  on-chip  electrodes 
to  the  nanobeam  is  of  utmost  importance  for  efficient  transduction  (Figure  1(c)  and  (d)).  Once 
free-standing  diamond  nanobeams  were  fabricated,  electrodes  were  patterned  on  the  diamond  sub¬ 
strate  via  a  standard  metal  lift-off  process.  First,  the  diamond  substrate  was  spin  coated  with  a 
polymethylmethacrylate-copolymer  (MMA/PMMA)  bilayer  resist,  where  the  MMA  copolymer 
thickness  was  carefully  chosen  to  be  slightly  thicker  than  the  distance  between  the  nanobeam  top 
surface  and  the  substrate.  Conformal  resist  coating  was  observed,  without  any  complication  due 
to  existing  device  geometry.  After  resist  coating,  exposure  and  alignment  were  done  via  electron 
beam  lithography.  50  nm  of  titanium  and  200  nm  of  gold  layer  were  then  evaporated  and  followed 
by  lift-off  in  remover  PG  to  complete  electrode  patterning.  Figure  2(b)  is  a  top-down  view  of  dia- 
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mond  nanobeam  cantilever.  It  can  be  seen  that  good  alignment  can  be  achieved,  with  error  on  the 
order  of  tens  of  nanometer  -  this  is  more  than  adequate  for  efficient  actuation.  In  fact,  we  note 
that  the  slight  misalignment  is  actually  beneficial  since  it  enables  the  actuation  of  in-plane  motion, 
as  discussed  below.  Figure  2(c)  shows  an  array  of  fabricated  diamond  doubly  clamped  nanobeam 
mechanical  resonators  that  share  driving  electrodes.  This  configuration  allows  us  to  characterize 
in  parallel  large  number  of  resonators  having  slightly  different  geometry  and  hence  different  reso¬ 
nance  frequencies.  Our  nanomechanical  resonators  had  wide  range  of  widths  (200  nm  -  300  nm) 
and  lengths  (1  [im  -  20  jum),  corresponding  to  fundamental  flexural  resonance  frequencies  (from 
Euler-Bernoulli  beam  theory)  ranging  from  a  few  MHz  to  hundreds  of  MHz.  Due  to  the  nature  of 
angled-etching  fabrication  technique,  width  and  thickness  of  nanobeam’s  triangular  cross-section 
are  correlated.  The  smallest  cantilever  that  could  be  characterized  (Figure  1(a))  had  a  width  of 
300  nm  and  the  length  of  4  [im,  while  the  smallest  doubly  clamped  nanobeam  (Figure  1(b))  had  a 
width  of  200  nm  and  the  length  of  7  gm. 

All  experiments  were  performed  at  room  temperature,  inside  a  vacuum  chamber  at  the  pressure 
below  10-5  Torr.  Figure  3(a)  shows  a  schematic  of  the  experimental  setup,  which  is  a  slightly 
modified  version  of  that  used  previously  for  characterizing  the  Brownian  motion  of  single-crystal 
diamond  nanobeams. 17  Specifically,  network  analyzer  was  used  instead  of  spectrum  analyzer  in 
order  to  obtain  driven  response  of  our  devices  (by  sweeping  the  driving  frequency  around  me¬ 
chanical  resonance),  and  a  bias-tee  was  included  to  combine  DC  bias  with  RF  drive  signal  to 
ensure  proper  actuation:  since  actuation  force  is  proportional  to  the  square  of  applied  voltage, 
F  oc  ( V Dc  +  Vrf  cos  cot)2. 19  DC  bias  is  needed  to  provide  force  component  at  the  driving  frequency. 

For  the  most  of  fabricated  diamond  nanobeams,  both  the  fundamental  out-of-plane  and  in-plane 
modes  were  characterized.  Resonant  responses  of  fundamental  out-of-plane  motion  of  devices 
shown  in  Figure  1(a)  and  (b)  are  plotted  in  Figure  3  (b)  and  (c),  respectively.  Solid  curves  are 
fitted  to  the  raw  data,  with  both  figures  showing  the  expected  resonant  responses  in  the  linear 
regime.  The  maximum  resonant  frequency  of  the  out-of-plane  mode  that  we  could  measure  using 
our  approach  was  ~40  MHz,  in  the  case  of  7  /im  long  doubly  clamped  nanobeam.  The  maximum 
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resonant  frequency  of  the  in-plane  mode,  measured  in  the  same  devices,  was  ^50  MHz.  To  the  best 
of  our  knowledge,  this  is  the  highest  actuation  frequency  of  flexural  mechanical  vibration  achieved 
by  dielectrophoretic  actuation  to  date.  Unfortunately,  in  our  current  experiments,  we  were  not  able 
to  measure  devices  with  resonances  >50  MHz,  due  to  the  deteriorated  signal-to-noise  ratio  of  our 
measurements.  Noise  floor  of  our  experiment  was  affected  by  three  different  instruments:  shot 
noise  from  laser  source,  dark  current  and  thermal  noise  from  photodetector  and  thermal  noise  from 
the  receiver,  which  can  be  either  the  network  analyzer  or  spectrum  analyzer.  Depending  on  the 
settings  of  instruments,  any  of  these  three  could  be  the  limiting  factor  for  the  noise  floor.  Most  of 
our  measurements,  however,  were  limited  by  the  shot  noise  from  the  laser  source. 

By  further  increasing  the  input  RF  power,  and  by  using  longer  devices  having  smaller  spring 
constant,  the  diamond  nanomechanical  structures  could  be  driven  deep  into  the  nonlinear  regime 
(Figure  4).  We  note  that  longer  devices  are  “softer”  and  therefore  deflect  more  for  the  same  applied 
voltage,  thus  easily  entering  nonlinear  regime.  Typically,  nonlinearity  in  small  scale  mechanical 
resonators  can  be  phenomenologically  modeled  using  Duffing  equation: 

where  x(t),£l o,  Q,  F(t )  and  m  are  the  beam  displacement,  the  resonance  frequency,  mechanical 
quality  factor,  external  driving  force  and  effective  mass  of  the  resonator,  respectively.  The  cubic 
term  f5x3(t)  in  the  equation,  so  called  ‘Duffing  nonlinearity’,  determines  the  nonlinear  behavior 
of  the  resonator:  when  /3  <  0  resonance  tends  to  move  to  the  lower  frequency  when  amplitude 
increases,  resulting  in  “softening”,  that  is  lowering  of  the  spring  constant;  the  opposite  is  true  when 
j8  >  0  and  spring  “hardening”  occurs.  Figure  4(a)  and  (b)  shows  softening  of  the  fundamental  out- 
of-plane  and  in-plane  modes,  respectively,  for  8  [im  long  diamond  nanobeam  cantilever.  On  the 
other  hand,  spring  hardening  is  observed  in  the  case  of  13  jim  long  doubly  clamped  diamond 
nanobeam,  for  both  out-of-plane  and  in-plane  modes  shown  in  Figure  4(c)  and  (d),  respectively. 

There  are  many  physical  origins  that  can  give  rise  to  mechanical  nonlinearities,  including,  but 


x(t)  = 


F(t) 


m 


(i) 
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are  not  limited  to  transduction  effects,  actuation  scheme,  material  properties,  non-ideal  boundary 
conditions,  damping  mechanics  and  geometric/inertial  effects.25  Based  on  our  FEM  simulations, 
nonlinearity  due  to  dielectrophoretic  actuation  can  be  ruled  out  due  to  small  deflection  (on  the  order 
of  nanometer)  of  the  resonator.  Nonlinearity  arising  from  transduction  on  nonlinearity  can  also  be 
ruled  out  due  to  the  relative  small  deflection  compared  to  the  laser  wavelength  used.  Geometrical 
nonlinearity  is  the  most  common  nonlinearity  observed  in  NEMS,  and  according  to  theory  results 
in  a  hardening  effect. 26  In  contrast,  inertial  nonlinearity,  another  physical  origin  of  nonlinearity, 
always  results  in  a  softening  behavior.27  Both  of  these  effects  can  be  significant  in  our  devices, 
and  their  interplay  determines  the  overall  sign  of  the  Duffing  constant  and  overall  softening  or 
hardening.26 

We  characterized  two  high  aspect  ratio  (length/width  >  20)  cantilevers  and  found  that  they  both 
showed  spring  softening  behavior  for  their  first  order  out-of-plane  and  in-plane  resonances.  Based 
on  these  results,  we  believe  that  the  nonlinear  behavior  observed  in  diamond  cantilevers  is  likely 
due  to  inertial  effects.  We  also  studied  ten  high  aspect-ratio  doubly  clamped  diamond  nanobeams, 
and  we  observed  that  the  sign  of  the  Duffing  constants  was  not  always  the  same.  While  in  all  cases 
fundamental  out-of-plane  modes  of  doubly  clamped  nanobeams  showed  spring  hardening,  the  fun¬ 
damental  in-plane  modes,  however,  showed  both  hardening  and  softening,  without  following  any 
particular  trend.  In  fact,  even  beams  of  nearly  the  same  geometry  on  the  same  chip,  revealed  oppo¬ 
site  mechanical  nonlinearites.  We  believe  that  theses  complexities  are  caused  by  residual  stress, 17 
that  was  apparent  in  some  high  aspect  ratio  beams  that  were  buckled  up  or  down  in  post-fabrication 
SEM  images.  The  sign  of  the  internal  stress  can  cause  different  types  of  nonlinearity:  for  compres¬ 
sive  stress  spring  softening  occurs,  while  for  tensile  stress  spring  hardening  occurs.28 

Nonlinear  behavior  was  further  confirmed  by  hysteresis  curves  as  shown  in  Figure  5(a)  and 
(b),  which  are  measured  from  fundamental  out-of-plane  and  in-plane  modes  of  19  ^tm  long  doubly 
clamped  diamond  nanobeam,  respectively.  This  device  was  particularly  chosen  among  others  as 
a  representative  examples  because  of  its  prominent  nonlinearity,  though,  other  devices  gave  the 
hysteresis  curves  as  expected.  Further  studies,  outside  the  scope  of  the  current  manuscript,  on 
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large  number  of  diamond  cantilevers  are  needed,  however,  to  unambiguously  identify  the  physical 
origin  of  nonlinearity  in  diamond  cantilevers. 

Our  dielectrophoretic  actuation  scheme  introduces  an  additional  effective  spring  since  the  actu¬ 
ation  force  has  dependence  on  position. 19  Furthermore,  since  the  force  has  a  quadratic  dependence 
on  applied  voltage,  the  mechanical  eigenfrequency  is  tunable  with  DC  bias,  and  the  amount  of  shift 
in  frequency  has  quadratic  dependence  on  voltage.  Figure  6  (a)  shows  Brownian  motion  measured 
with  a  spectrum  analyzer  as  the  applied  DC  bias  was  changed  from  —9  V  to  +9  V.  Bright  spots  of 
each  data  column  correspond  to  the  Brownian  motion  resonance  location.  The  solid  black  line  is 
a  quadratic  fit  for  applied  DC  bias  and  shows  an  excellent  match  with  the  theoretical  prediction. 
In  the  given  range  of  applied  DC  bias,  the  mechanical  resonance  could  be  tuned  over  roughly  260 
full  width  at  half  maximum  of  the  zero  bias  Brownian  motion  peak.  We  observe  a  blue  shift  of 
resonance  frequency  because  this  motion  is  in-plane  vibration  and  the  effective  spring  has  always 
the  same  sign  with  the  elastic  spring.  In  the  case  of  out-of-plane  modes,  frequency  tuning  can  be 
either  red  or  blue  shifted,  depending  on  the  height  of  the  beam  from  the  driving  electrodes.20 

Since  the  resonance  frequency  is  easily  parametrically  tuned  more  than  twice  the  linewidth, 
parametric  excitation  is  also  expected.  Criteria  for  a  resonator  to  be  used  for  parametric  actuation 
is  given  in  equation  (2). 29 


d£ 

dV 


■V  ■  —  >  2 

v  max  r  ^  ^ 
JO 


(2) 


where  /,  V  and  Q  are  resonance  frequency,  amplitude  of  applied  RF  voltage  and  mechanical  quality 
factor,  respectively,  /o  is  a  resonance  frequency  without  RF  voltage  input  and  Vmax  represents  the 
maximum  voltage  that  can  be  applied.  Since  the  spring  constant  is  a  function  of  position,  it  can 
be  modeled  with  Mathieu’s  equation  with  a  Duffing  nonlinearity.  Equation  (1)  may  be  modified  to 
include  parametric  tuning  as  shown  below: 


a2 


—jp  +  T  (l  T  oc  +  /3x2(t)  —  2T sin2£2ot) 


x(t)  =  0 


(3) 


where  a  is  the  detuning  between  the  parametric  excitation  and  the  fundamental  mode  £2o,  T 
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is  proportional  to  the  parametric  excitation  amplitude,30  whereas  other  parameters  are  defined  the 
same  way  as  in  equation  (1).  Right  hand  side  of  the  equation  (3)  is  zero,  since  there  was  no  direct 
driving  at  resonance  frequency.  When  the  resonator  with  eigenfrequency  /o  can  be  parametrically 
excited,  its  excitation  frequency  can  be  for  any  integer  n.  In  most  of  the  cases,  response 
at  «  =  1  is  used  since  it  is  the  strongest,  although  the  submultiples  of  it  has  been  observed  in 
experiment.31  Mathieu’s  equation  with  Duffing  nonlinearity  (3)  can  be  analytically  solved  and  the 
solution  predicts  its  stability  on  a  phase  plane,  axes  of  which  are  detuning  and  driving  amplitude. 
Here,  we  show  “instability  tongue”32  when  doubly  clamped  diamond  nanobeam  is  parametrically 
excited  with  n  =  1.  In  Figure  6  (b)  the  measured  instability  tongue  is  shown  when  nanobeam 
was  excited  around  the  twice  of  its  natural  frequency  of  7.25  MHz.  In  this  experiment,  excitation 
was  done  by  RF  signal  generator  and  the  response  was  measured  with  spectrum  analyzer,  with 
the  amplitude  of  driven  motion  was  calculated  by  simple  signal  processing.  Parametric  excitation 
is  particularly  interesting  for  NEMS  devices  since  it  can  circumvent  electric  cross  talk,  which 
can  be  detrimental  for  nanoscale  system,33  and  can  be  used  to  realize  a  NEMS  oscillator29  and 
mechanical  memory  element. 30  Furthermore,  parametric  oscillators  are  of  interest  for  mechanical 
quality  factor  enhancement,34  parametric  amplification,35  and  noise  squeezing.26 

In  summary,  we  have  realized  a  resonant  actuator  based  on  dielectrophoresis  for  single-crystal 
diamond  nanomechanical  systems,  thus  successfully  demonstrating  single-crystal  diamond  NEMS. 
Actuation  of  both  cantilever  and  doubly  clamped  diamond  nanobeams  was  achieved  for  both  the 
fundamental  out-of-plane  and  in-plane  vibrations.  Driving  frequency  range  spans  from  a  few  MHz 
to  nearly  50  MHz,  though  higher  frequency  actuation  is  expected  to  be  realized  by  improving 
the  signal-to-noise  ratio  of  the  optical  interferometric  displacement  detection,  or  moving  to  other 
low-noise  displacement  read-out  schemes.  The  dielectrophoretic  actuation  scheme  employed  is 
expected  to  maintain  the  mechanical  quality  factor  of  the  base  diamond  device,  and  also  to  not 
compromise  the  optical  properties  of  integrated  diamond  color  centers. 

Furthermore,  our  actuation  scheme  can  efficiently  drive  the  mechanical  resonance  well  into  the 
nonlinear  regime,  with  reasonable  RF  power  levels.  Studying  physical  origins  of  mechanical  non- 
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linearity  in  NEMS  devices  are  active  field  of  research,  still  in  its  infancy  with  many  open  questions 
unanswered. 36-38  To  the  best  of  our  knowledge,  we  believe  that  our  work  is  the  first  demonstra¬ 
tion  of  nonlinear  mechanical  response  of  single-crystal  diamond  resonator.  Single-crystal  diamond 
has  a  great  potential  as  a  nonlinear  NEMS  platform  due  to  its  ultrahigh  mechanical  quality  factor. 
Recently,  it  has  been  shown  that  single-crystal  diamond  resonators  can  have  mechanical  quality 
factor  over  300,000 18  or  even  over  a  million16  at  room  temperature.  In  this  regard,  the  capability 
of  our  system  to  preserve  high  mechanical  quality  factor  is  of  great  importance  in  the  study  and 
application  of  nonlinear  mechanics.  Although  pre-stressed  nitrides  in  general  show  the  highest 
quality  factor  in  NEMS,38  it  is  not  suitable  to  study  nonlinear  mechanics  because  of  its  strong 
tensile  stress.26 

Lastly,  our  actuation  scheme  is  capable  of  strong  frequency  tuning  simply  by  applying  DC  bias, 
and  leads  to  useful  parametric  excitation.  Parametric  excitation  is  expected  to  be  important  in  many 
applications  including  NV-  center  engineering,  since  actuation  signals  can  be  completely  filtered 
out  from  desired  NV-  center  manipulation/read-out  signals  in  frequency  domain.  Single-crystal 
diamond  actuation  we  developed  here  is  expected  to  be  a  gateway  to  a  more  sophisticated  platform 
for  the  control  of  NV-  center  via  nanoscale  mechanics. 
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Figure  1:  SEM  images  of  (a)  4  jiim  cantilever  and  (b)  7  /tin  doubly  clamped  beam.  These  rep¬ 
resent  the  shortest  devices  that  could  be  characterized  to  date,  (c)  Finite  element  method  (FEM) 
simulation  is  used  to  calculate  applied  force  to  suspended  nanobeam  with  a  given  geometry  and 
electrostatic  environment.  Color  map  stands  for  the  voltage  and  the  streamlines  show  correspond¬ 
ing  electric  field,  (d)  Vertical  force  per  unit  length  applied  to  the  beam  from  (c)  is  plotted  as  a 
function  of  beam  width  and  height,  in  which  applied  voltage  of  20  V  was  assumed.  It  can  be  seen 
that  force  is  maximized  when  height  is  small.  In  fact,  there  is  an  optimal  position  that  mazimizes 
the  force,  since  force  goes  to  zero  in  the  middle  of  two  electrodes  due  to  the  symmetry  and  it  also 
goes  to  zero  at  the  far  distance  due  to  the  lack  of  electrical  field  gradient. 


16 


Figure  2:  (a)  Schematic  illustration  of  angled-etching  nanofabrication  of  for  suspended  diamond 
nanobeams  and  the  following  standard  bi-layer  MMA/PMMA  process  for  lift-off.  Each  step  is 
composed  as  the  following:  (i)  Electron  beam  lithography  mask  is  deposited,  (ii)  top-down  reac¬ 
tive  ion  etching  of  diamond  is  performed,  followed  by  the  (iii)  angled-etching  step  and  (iv)  mask 
removal,  (v)  New  e-beam  resist  is  deposited,  and  (vi)  electron  beam  lithography  followed  by  (vii) 
metal  evaporation  and  (viii)  lift-off  are  used  to  define  electrodes,  (b)  High  mag  SEM  image  of 
7  jum  cantilever  shows  that  good  alignment  can  be  achieved,  on  the  order  of  tens  of  nanometers, 
which  allows  for  efficient  actuation,  (c)  SEM  image  of  device  array  sharing  electrodes. 
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Figure  3:  (a)  Optical  characterization  setup.  Fundamental  out-of-plane  resonant  response  of  de¬ 
vices  shown  in  Figure  1(a)  and  (b)  are  given  in  (b)  and  (c),  respectively. 
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Figure  4:  When  nanobeams  are  driven  hard  enough,  they  start  to  show  nonlinear  behavior.  Nonlin¬ 
ear  response  of  8  [im  long  cantilever  for  its  (a)  out-of-plane  and  (b)  in-plane  modes.  Softening  of 
spring  was  observed  for  both  types  of  modes.  On  the  other  hand,  spring  hardening  was  observed 
in  the  case  of  strongly  driven  13  [im  long  doubly  clamped  diamond  nanobeam  both  for  its  (c) 
out-of-plane  and  (d)  in-plane  modes. 
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Figure  5:  Hysteresis  curve  of  (a)  hardening  and  (b)  softening  mechanical  nonlinearity  measured 
from  19  jUm  long  doubly  clamped  diamond  nanobeam  when  frequency  was  swept  up  and  down, 
(a)  and  (b)  are  the  fundamental  out-of-plane  and  in-plane  modes,  respectively. 
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Figure  6:  (a)  Tuning  of  mechanical  resonance  using  DC  bias.  With  applying  ±9V,  frequency 
tuning  range  that  can  be  achieved  is  approximately  260  linewidths.  Quadratic  nature  of  tuning  can 
be  used  to  control  the  onset  of  parametric  instability  in  terms  of  RF  power,  (b)  Typical  tongue 
shape  of  parametric  instability  was  observed.  Low  response  spots  inside  the  instability  region  are 
attributed  to  data  acquisition  software  problems. 
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